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^\ ' We study the effect of disorder and interactions on a recently proposed magnetic field induced 

. localization mechanism. We show that both partially destroy the extreme confinement of the exci- 

tations occuring in the pure case and give rise to unusual behavior. We also point out the role of 
the edge states that allows for a propagation of the electrons in these systems. 
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I. INTRODUCTION 



The dynamics of a quantum particle in a periodic potential in the presence of a uniform static magnetic field has 
revealed many beautiful effects and has been the subject of ongoing researches for several decades. The competition 
' between the two characteristic length scales involved, namely the lattice period and the magnetic length produces a 
. very conuplex pattern of energy levels. In twij dimensions, this leads to fractal structures in the spectrum anticipated 
O by AzbelEl, and further studied by HofstadterEl, Claro and Wanniercl, Rammalo, and many others. Some striking experi- 
, O , . mental manifestations of these effects have been found in macroscopic properties of two-dimensional superconducting 
wire networks, such as the variation of the critical temperature with respect to the external magnetic field, the 
magnetization or the critical current. The connection between the linearized Ginzburg-Landau equations used to 
describe these wire networks close to the superconducting transition, and the tight-binding spectrum of electrons had 
been established slightly before these experiments by De GewiesB and AlexanderQ. 

More recently, a surprising effect has been presented in refJ3 in a two-dimensional lattice with hexagonal symmetry, 
the so-called Ta lattice displayed in Fig. For some values of the magnetic field corresponding to half a flux quantum 
per elementary tile, the energy spectrum of a tight-binding model with nearest neighbor hopping collapses into three 
highly degenerate levels. Furthermore, we have shown that it is possible to build energy eigenstates where the 
^ , probability to find an electron is non vanishing only in a finite size cluster which we have called an Aharonov-Bohm 
^ ' cage. This corresponds to a localization mechanism due to quantum interferences of Aharonov-Bohm type between 
paths enclosing a half-integer number of flux quanta. For instance, if an electron is initially located on a given site 
^ : of the lattice, it will never propagate beyond the boundary of the cage associated to this site, provided the magnetic 
' field is chosen properly, as discussed above. 

Since this theoretical observation, several experiments have searched for manifestations of these cages, on real 
systems. Aj-serie of investigations on superconducting wire networks has been completed by Abilio et al.u and by 
^ , Serret et ala providing a good agreement between predictions and measurements of the critical temperature and the 
critical current. In addition, the unusual nature of the mixed state for half a flux quantum per elementary tile has 
been clearly indicated by a strong reduction ftfi-^he critical current, and by magnetic decoration experiments which 
shows a very highly disordered vortex patternO'tj. 
5— i A slightly more direct probe of this localization effect has been provided by transport measurements in, mesoscopic 

lattices carved in a two-dimensional electron gas at the interface of GaAs/GaAlAs heterostructuresEJ. In theses 
systems, the number of tranverse conduction channels in the wires can be very small (a few units), and the mean 
free path much larger than the distance between two nearest neighbor nodes. At intermediate magnetic fields, Naud 
et al. have clearly observed a periodic modulation of the magnetoconductance with a period of one flux quantum 
per plaquette (h/e). In the same experimental conditions, a square lattice does not show these oscillations. It seems 
more than plausible that these /i/e-periodic oscillations seen on lattice could be attributed to the Aharonov-Bohm 
cages. Indeed, this localization mechanism establishes a strong difference in the tranport properties between integer 
and half integer fluxes, which is expected to survive in the presence of a weak enough disorder as discussed in ref.E3. 
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These experiments clearly motivate us to address the question of the robustness of these cages in the presence 
of "real life" perturbations such as disorder, finite size effects and electron-electron interactions. The main purpose 
of this paper is to attempt to bridge the gap between the very idealized tight-binding model initially studied in 
ref.l3, and the experimental systems mentioned above. As a general trend, we will show that if the cages are in 
themselves very fragile, their existence in the ideal system induces a few remarkable features for the perturbed ones. 
The underlying reason for this is that these perturbations act on a very degenerate system, and as often in physics, 
degenerate perturbations are expected to induce some very interesting phenomena. One of the most famous example 
is, of course, the fractional quantum Hall effect where the degeneracy of the kinetic energy at fractional fillias; factors 
opens the way to the formation of strongly correlated many-body exotic states such as the Laughlin liquidO. 

This paper is organized as follows. Section || presents a detailed discussion of the butterfly-like energy spectrum 
for a tight-binding model on the T3 lattice. Section [II focusses on half integer fluxes for which the cages appears. 
Several viewpoints are given to provide a simple physical intuition of what happens in the system. This enables us 
to show other examples of lattices for which Aharonov-Bohm cages exist. Section IV discusses disorder, effects still 
within the tight-binding framework. So, it is complementary to the previous study presented in rcf.EJ where the 
continuous wire model was analyzed, in connection to the experiments on two-dimensional gases. Here, we investigate 
some properties of one particle eigenstates such as the inverse participation ratio as a function of the disorder strength 
and the magnetic field. Section ^ is devoted to edge states. These states are found to be confined in a finite width 
strip along the boundaries of the sample. We outline the subtleties involved in the translation from the tight-binding 
model properties to those of the continuous one-dimensional wire networks. Finally, section Vl adresses interactions 
effects in the context of the Hubbard raodel. Most of the results deal with the two electron case. In agreement with 
a previous study for the chain of loopscfl, we find that interactions are able to create some extended and dispersive 
two particle eigenstates. This can be viewed as a kind ofjjiteraction induced delocalization phenomenon, similar to 
the one discussed by Shepelyansky for disordered systemsta. The presence of the cages in the non interacting case is 
reflected by the fact that these extended states involve a close binding of the two particles in real space, even if the 
local interaction is repulsive. We also give some ground state configurations for a finite density of particles and we 
show that the energy per particle has a singular behaviour when the number of electrons per site becomes larger than 
1/3. In the last section, we conclude and discuss the experimental relevance of our work. Technical details for the 
two interacting particles problem can be found in the appendix 



II. SPECTRAL PROPERTIES 



We consider the 73 lattice displayed in Fig. |^ which is a bipartite periodic structure with 3 sites per unit cell: one 
6-fold coordinated site A and two 3-fold coordinated site B and C. 




FIG. 1. A piece of the Ts lattice embedded in the perpendicular magnetic field H ; a denotes the lattice spacing. 
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This tiling can also be seen as the dual of the most famous Kagome lattical3. We consider a tight-binding hamil- 
tonian defined by : 

<jj> 

where |i) is a localized orbital on site i. When_H = 0, the hopping term = 1 if i and j are nearest neighbors and 
otherwise. In the presence of a magnetic fieldEII, tij is multiplied by a phase factor e'^*^ involving the vector potential 
A : 

7„ = — / A.dl, (2) 

where = hc/e is the flux quantum. In the following, we focus on the case of a uniform magnetic field H = Ht. 
which can be obtained, for example, with the Landau gauge A = 7J(0, a;,0), and we denote by </> = Ba?\/?>/2 the 
magnetic flux through an elementary rhombus. 

Let us first briefly analyze the zero field case H = 0. In this case, the one-particle spectrum is given by the following 
dispersion relation : 

E± (k) = ± ^6 + 4 [cos(k.u) + cos(k.v) + cos(k.w)], (3) 

where u = a(3/2, -\/3/2) and v = a(3/2, — \/3/2) are the vectors of the primitive cell (see Fig. |l|), w = v — u, and 
where k is a wave vector lying in the first Brillouin zone. The associated eigenstates are standard Bloch waves. In 
addition, there is a non dispersive band at Eq = resulting from the bipartite character of the structure, and whose 
degeneracy is given by the difference between the number of 3-fold and 6-fold coordinated sites. Since its origin is 
purely topological, this energy will always_be an eigenvalue for any H, with the same degeneracy. The density of 
states, originally computed by Sutherlandts, is shown in Fig. 
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FIG. 2. Density of states of the Ts lattice for H = 0. 



For H ^ 0, the spectrum is of course much more complex since there remains only one translational invariance 
along the y direction due to the gauge choice, so that the eigenfunctions can be written as : 



(Pj{x,y) ^ 'ijj.j{x)e^''^y, ky e 0,27r/a\/3 
where the index j = A, B ,C . Thus, the secular equations read : 
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where k = and 7 = 27r(/)/0o = 27r/. Note that with the origin O chosen in Fig. x only takes integer or 

half-integer values. If one looks for solutions e ^ of the secular system, one can substitute and (|^) in (||) to 
obtain an effective one-dimensional equation that only involves six-fold coordinated sites : 



(e^ - 6)?A™ = 2 cos <^2cos 
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where ipm — ipAix = 3ma/2), (to e Z). As it can be readily seen in Fig. the six-fold coordinated sites form 
a triangular lattice so that Eq. (||) has to be compared to those derived by Claro and Wannier for the triangular 
latticed : 
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with similar notations. So, the lattice spectrum can be simply obtained from the triangular lattice spectrum by 
setting 7t = 37/2, since one has : 



- 6 = 2 cos 
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ST- 



(10) 



Since the tiling is bipartite, the spectrum is symmetric (e <-*■ —s). Moreover, Eq. (^) displays a translation symmetry 
/ /+ri(Vn € Z) and a reflection invariance about half-integer values of/. We thus limit our analysis to < / < 1/2. 
For rational values oi f = p/q {p,q mutually prime), the system (^) becomes closed after a translation by q periods 
and the spectrum is made up of 2q bands. Note that, when q — 3q' {q' € N), this period is actually reduced by a 
factor 3 and the spectrum has only 2q' bands. In particular, for / = 1/3, one can exactly compute the dispersion 
relation : 

£^(k, fc^;) — 6 = — 2 cos(2k) + 4 cos(k) cos(A:2;a) , (11) 

where (fc^, k) € [0, 27r/a\/3] x [0, n]. For this flux, the spectrum spreads from —3 to -1-3 and is gapless. Figure || shows 
the 7^ lattice spectrum support versus the reduced flux /. The most spectacular and unusual feature is that, for 
/ — 1/2, the spectrum collapses into three eigenvalues £0 = and e± — ±V6 as it can be readily seen from Eq. (px|). 



(8) 




FIG. 3. Spectrum of the T3 lattice as a function of the reduced flux. 



This is all the more curious that usually, for an infinite periodic structure and for rational values of /, the spectrum 
is absolutely continuous (band- like). Here, it behaves as a super-atom with three infinitely degenerate levels with 
equal spectral weights 1/3. The existence of these non dispersive bands suggests the possibility to build Wannicr- 
type localized eigenstates. This remarkable fact is, for any electron density, susceptible to induce original physical 
properties reminiscent of those of a localized system. 
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III. THE AHARONOV-BOHM CAGES 



In this section, we shall analyze these properties from the quantum dynamics point of view by characterizing 
the spreading of a wave packet in the Ts lattice, at / = 1/2. To achieve this, it is worth focusing on the spectral 
charateristics at a local level. The magnetic field being uniform, the whole spectrum is indeed recovered from the local 
density of states (LDOS) on the two different types of sites. At this point, it is more convenient to shift from the above 
Landau gauge to a cylindrical (symmetric) gauge defined by A = H{—y/2, x/2, 0) that respects the discrete rotational 
symmetry existing locally. We proceed to an analytic Lanczos tridiagonalization using the recursion methodt^. The 
principle of this algorithm is to generate an semi-infinite chain according to the recursive relation : 



which allows to evaluate the LDOS on the initial orbital |<^o)- The diagonal and off-diagonal terms (a„,6„) are 
obtained by requiring the normalization of the orbitals \^n)- Note that in our case, all the diagonal elements a„ 
vanish because : (?) the structure is bipartite, (ii) Hq is purely off-diagonal. 

Let us first compute the recursion chain associated to a six-fold coordinated site and choose, as initial orbital, 



The first recursion orbital \tpi) is a linear symmetric combination of the first shell sites (three- fold coordinated grey 
sites). At the next step, one obtains : 



bn+l\fn+l) = {Ho - a„)|(^„) - bn\fn-l 



(12) 



|(po) = |0) (see Fig .|). One has : 




(13) 




(14) 
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FIG. 4. Aharonov-Bohm cage associated to the six- fold coordinated sites and its recursion chain. For / = 1/2, the black 
sites are never visited by a wave packet initially localized on the central site (white). 

For / = 1/2, 62 = so that the LDOS is reduced to two values e± = ±\/6 with an equal spectral weight 1/2. This 
effect can be simply understood in terms of Aharonov-Bohm effect. Indeed, the amplitude of probability ►2-»9 to 
go for example, "in two steps", from to 9 via 2, is exactly the opposite of ^o^3^9j so that the resulting amplitude is 
zero. Then, any wave packet initially localized on a,six-fold coordinated site = 0)) = |0)) is completely trapped 
inside what we have called an Aharonov-Bohm cageQ whose precise definition is given below. One can easily determine 
its time evolution that is periodic and given by : 

!■(/'(*)) = cos 

(/3o) ~ i sin (Vet) . (15) 

In particular, the autocorrelation function |(0|i/'(t))p = cos^ (\/6t) for / — 1/2, whereas for / = : 



(16) 



{IBZ stands for first Brillouin zone) and behaves at large t as : 
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We emphasize that for a generic rational / 7^ 1/2, the hamiltonian is periodic and one expects |(0|'(/'(t))P ~ l/i^- 
A contrario, when / is irrationnal, the hamiltonian is quasiperiodic and one expects, as for the square latticeE3, 
|(0|'(/'(t))P ~ l/t"^" with a < 1 since the spectrum should be singular continuous. 

The same type of recursion calculation can be made for the three-fold coordinated sites. In this case, one easily 
obtains : 



61 = \/3, 



62 = v/4cos2(7r/)+3. 



(18) 
(19) 
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8cos2(7r/) + 4 [cos(^/) + cos(37r/)]" 
^^-y 4cos2(vr/)+3 ' 

so that 63 = for / = 1/2, and the LDOS on this type of site reduces to the £0 = (spectral weight 1/2), and 
e± = ±-\/6 (spectral weight 1/4). The quantum evolution of a wave packet initially localized on such a site is 
therefore, as previously, confined in a cage displayed in Fig. [s] that is larger than the one obtained for the 6-fold 
coordinated sites. 




FIG. 5. Aharonov-Bohm cage associated to the three-fold coordinated sites and its recursion chain. 



Finally, since the Ta lattice is only made up of 6-fold and 3-fold coordinated sites, any wave packet (with a finite 
initial extension) will be confined inside a bound Aharonov-Bohm cage at / = 1/2. We emphasize that this localiza- 
tion phenomenon must be understood from the quantum dynamics point of view. Indeed, contrary to the Anderson 
localization, the eigenstates of the system are not exponentially localized, and not even localized at all since the high 
degeneracy allows one to build any type of eigenstates (possibly extended). 

To proceed further, it is necessary to give a precise definition of these Aharonov-Bohm cages. For a given wave 
packet submitted to a magnetic field characterized by a reduced flux /, we define the Aharonov-Bohm cage 

as the set of sites visited by this initial wave packet during its evolution. In general, is infinite, but as shown 
previously, it can, for specific values of the magnetic field (/ — 1/2 for the lattice), be bound. To analyze the 
cage structures for any independent electron model, it is sufficient to characterize the cages of all inequivalent sites. 
Indeed, the superposition principle implies that if \^p{t = 0)) = X]i=i '^i '^^^ following property : 

n 

C U C/ . (21) 

i=l 

With this definition, several cases can be encountered : 

• all cages are unbound for any / (ex : square lattice, triangular lattice, honeycomb) ; 

• some cages are bound for particular values of / (ex : the Penrose tiling and the octagonal tiling displayed in 
Fig. I) ; 

• all the cages are bound for the same values of / (ex : T3 lattice). 
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FIG. 6. Two examples of structure presenting bound cages for special initial conditions of the wave packet localized at the 
center of some 5- fold symmetric stars for the Penrose tiling (left) and 8- fold symmetric stars for the octogonal tiling (right). 

This latter case, that we will call fully confined structure, can be met for other tilings. As an example, we have 
displayed in Fig. ^ , the so-called lattice which can be obtained by deforming an approximant with 14 sites per 
unit cell of the octagonal tiling. 




FIG. 7. A piece of the periodic T4 tiling before (left) and after deformation (right). The five different types of sites are 
represented {»,o)- 

This structure has exactly the same physical properties as the % lattice at / = 1/2 (/ refers here to the reduced 
flux measured in the smallest tile). Thereafter, we shall consider a light version of this tiling obtained by eliminating 
the sites o and their related bonds (see Fig. 0). The resulting structure is simply a set of connected 8- fold symmetric 
stars with two different lengths whose ratio is equal to ^/2. We thus allow for two different hopping terms and ts 
for long and short edges respectively. For / = 1/2, it is always possible to choose a gauge such that all the tij's are 
real. A possible choice is represented in Fig. H. 
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FIG. 8. A possible gauge choice for the T4, lattice at / = 1/2. All the Peierls terms are equal to +1 except those represented 
by dotted lines that are equal to —1. For clarity we have omit the hopping term tL,s- 

Note that this gauge has the same periodicity as the structure. This is quite surprising since it is not the generic 
case (see next section for the 7^ lattice). It is then easy to diagonalize Hq and to obtain the thirteen following 
eigenvalues : (three- fold degenerate), ±^j2ij^ + 2ti. (two- fold degenerate), ± \/ts, + and ±^JAij^ + 2tg. 

For any hopping terms tL,s these eigenvalues do not depend on any wave vector and thus form the exact analogous 
of the non dispersive bands obtained in the Tjj lattice. 

We would like to point out that fully confined systems can also be obtained in one dimension. We have shown 
in Fig. ^ a quasi-one-dimensional structure that clearly displays bound Aharonov-Bohm cages at / = 1/2. Its 
eigenspectrum is, as for the % lattice, made up of three non dispersive bands with energies 0, ±2. 




FIG. 9. The simplest example of a periodic fully confined structure for / — 1/2. 

After this detailed description of this magnetic field induced localization, it is natural to wonder if this phenomenon 
is robust to various type of "perturbations" . For example, a non uniform magnetic field has been shown to drastically 
chanae. the energy spectrum and the localization properties of the system providing a complete destruction of the 
cagescil. In the next section, we study the effect of different kind of disorder and we show that, the destructive 
interference that leads to the Aharonov-Bohm cages are partially destroyed. 

IV. THE DISORDER 

A first natural way to introduce disorder in a system is to directly modify its structure by incorporating defects. 
In the case of the T3 lattice, such defects can be simply generated by a local elementary flip of three tiles as shown in 
Fig. 0. 




FIG. 10. Elementary flip for two-dimensional tilings. 

For / = 1/2, it is clear that a finite number of flips does not modify the bound nature of the cages since, sooner or 
later, a given wave packet will be embedded in the "pure" T3 lattice geometry which is responsible of the confinement. 
Nevertheless, one can wonder if this feature still holds for a finite density of such defects. Indeed, if we consider a 
large number of flips in order to get closer and closer to a random tiling configurationE3, we expect the destructive 
Aharonov-Bohm interferences to completely disappear. It would be, in this case, interesting to determine the critical 
density below which all the cages remains bound. This could be achieved by considering an extended initial state and 
by computing its time evolution for different realizations of the disorder and for different flip densities. Note that this 
value should be easily understood in terms of cages percolation threshold. 

Another possible way to introduce disorder is to randomly modulate either off-diagonal or diagonal terms in the 
hamiltonian. In this study, we mainly focus on the latter case and consider the following hamiltonian : 

i/d... = , (22) 
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where the tij are defined as in and where the on-site energies Vi are mutually independent gaussian random 
variable with variance a^. 

As the strength of disorder increases, the strongly localized cage eigenstates begin to spread over larger distances. 
In some sense, a weak disorder increases the one-particle localization typical length ! However, since the system is 
two-dimensional, we expect the overall localization length to be finite for any value of u 7^ and for any energy. For a 
given reduced flux /, as a is increased from zero, the one-particle eigenstates first get mixed within the energy bands 
of the pure model. This regime corresponds to the standard Anderson localization. This lasts until a becomes of the 
order of some energy gaps allowing then for interband mixing. 

For / = 1/2, the remarkable structure of the spectrum at zero disorder is responsible for the absence of the first 
disorder induced localization regime. Therefore, we expect, in this case, that most observables will depend only weakly 
on tJ, as long as a remains smaller than the gaps. To analyze this problem, it is interesting to construct an effective 
model in the lowest energy subspace of Hq, expressed in the localized cage basis (see Fig. |l^ in the next section). 
Denoting the normalized cage eigenstates associated to e_ = — \/6 centered around the six-fold coordinated sites r by 
|r), this effective hamiltonian has the form : 

if.//. = E^~'-.r'|r)(r'|+E^r|r)(r|, (23) 

(r,r') r 

where (r, r') stands for the sum over nearest neighbors in the triangular lattice generated by the the six-fold coordinated 
sites. The effective parameters fr,r' and Vr depend on the Vj^s. For simplicity, we use the indexation of Fig. ^so that : 

Vo = ~V6 + ^ + ^j2^. (24) 

io.8 - ^(Fi - V^2)e*"°'« , (25) 

where ao,8 is a gauge dependent phase factor. The general patterns for Vr and tr,r' can be easily inferred from 
these examples. We observe that Heff. involves both diagonal and off-diagonal disorder with comparable magnitude 
since : (V^'^) — — 7cr^/24 and (|ir,r'P) — (|ir,r'|)^ = ct^/72. It is interesting to note that non vanishing tunneling 
amplitude from one cage to a neighboring one are induced by an asymmetry between the two diagonal energies of the 
intermediate 3-fold coordinated sites connecting the cage centers. According to this effective model, the statistical 
properties of the eigenstates and the energy levels are independent of cr, provided the eigenenergies are rescaled 

We have checked this simple picture by numerical diagonalization of finite clusters of size as large as 605 sites 
with open boundary conditions. As discussed in the next section, this introduces edge states. Here, the shape of the 
clusters has been properly chosen to minimize the number of these eigenstates whose energy, in the pure system at 
/ = 1/2, differs from 0,±-\/6. To characterize the degree of localization of a given eigenstate 1-0), it is convenient to 
compute the inverse participation ratio (IPR) defined by : 

IPR(V') = ^'^'^'^^ . (26) 
(EJ^.P) 

This quantity reaches a non vanishing constant value in the large system size limit for localized states whereas it 
behaves as l/Ng for extended states, Ng being the total number of sites. For each value of a, we have averaged over 
10^ realizations of the disorder and computed the average value of IPR over the whole energy spectrum. The results 
are plotted in Fig. |ll| for several values of the reduced flux, as a function of the disorder strength a. 
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FIG. 11. Variation of the averaged log (IPR) as a function of the disorder strength a, for several values of / {Ns — 605 sites). 
The average is done over energies and over disorder realizations. 



We first notice that, for / = 1/2 the averaged IPR is indeed almost independent of cr up to cr ~ 1, in agreement 
with the qualitative picture given before. As a is further increased, the averaged IPR increases which is interpreted 
as a relative relocalization associated to an interband mixing. For / < 1/2, another energy scale D{f) related to the 
mean band width in absence of disorder emerges. If cr < D{f), the localization length is finite for any a and decreases 
with (T as indicated, for instance at / = 0.495, by the corresponding increase of the IPR for low a. This lasts until 
D{f) o" ^ 1 where eigenstates no longer evolves with a. The intermediate regime in smoothly connected to the 
plateau observed for f = 1/2. It is, in some sense, what remains of Aharonov-Bohm cages as the field is no longer 
fixed to its special value / = 1/2. Finally, the strong disorder regime (cr > 1), where the localization length becomes 
comparable to the lattice spacing, is found not to depend sensitively of / which sounds quite reasonable. 

Similar conclusions can be drawn by looking at the density of states. Plots of this quantity at / = 1/2 for various 
values of a are shown in Fig. ^ Note that the splitting of the three main peaks for a = 0.01 is due to edge states. 
The peaks broaden linearly with a in the plateau region which clearly ends when they merge. 

Finally, we would like to mcption a related study concerning the influence of the disorder on the transmission 
properties of quantum networkal3. For sufficiently weak disorder, it has been shown that the periodicity with respect 
to the magnetic flux of the magnetoresistance remains equals to h/e in the T3 lattice whereas it is completely dominated 
by the weak localization regime {h/2e) for more conventional structures, e. g. the square lattice. This indicates that 
the Aharonov-Bohm cages resists to a small amount of disorder. We emphasize that beautiful experiments by Naud 
et al. have recently shown up this phenomenon in two-dimensional electron gas (GaAs/GaAlAs)", measuring the 
magnetoresistance of mesoscopic artificial structures with the T3 geometry. 
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FIG. 12. Variation of the averaged density of states at / = 1/2 for various disorder strength a {Ns — 605 sites). 



V. THE EDGES STATES 



As already mentioned, interesting edge state properties appear a,t f — 1/2. We could study these states explicitely 
for the tight-binding hamiltonian Hq but, for the sake of comparing with the experimental situations, we shall consider 
a continuous model where the network is made up of perfectly conducting one-dimensional wires. The Schrodinger 
equation for a free particle of mass m moving along a wire reads : 

where x denotes the coordinate along the wire. For simplicity, we assume that no magnetic field is present but we will 
introduce it as soon as we will need it. If the extremities of the wire corresponds to nodes i and j, and if the origin 
is chosen such that x — for node i and x — a for node j, we may write the general solution along the wire as : 

sm(ka) 2m 

For clarity, we have omit the index fc, and we have set ipi — ip"^^ {x = 0) and ipj = ip"^^ (x = a). The amplitudes on 
the nodes are constrained by a linear boundary conditions at each node of the network. A frequently used condition 
which is compatible with current conservation at the node i is : 



E 



dip'^ix) 



= . (29) 



Here, the sum is taken over the Zi wires emerging from the node i. If we also impose the continuity condition of the 
wave function at each node {ip^^ {x = 0) independent of j), we obtain the following set of equations : 

where we have introduced A — cos ka. So, a single value of A generates an infinite discrete family of eigenstates with 
energy E{kn^x) — ^ {kx + provided A = cos kxa. To describe each eigenstate only once, we may impose k 

to lie in the interval [— 7r/a, tt/o]. Note that Eq. ( |30|) cannot be directly interpreted as an eigenvalue equation for 
a tight-binding problem since, in the general case, Zi varies from one site to another. For an infinite T3 lattice, or 
a finite one with appropriate boundary conditions, the coordination number Zi can only be equal to 3 or 6. It is 
therefore possible to perform the same decimation of the three fold-coordinated sites as for the tight-binding problem 
to obtain the same equation as Eq. (§) by setting — 18 A^. This simple correspondence between both problems 
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also holds in the presence of a uniform magnetic field. In zero field, since e runs from — \/T8 to vTS, A runs from 
— 1 to +1 and therefore k\ from [— 7r/a, 7r/a]. So, there is no gap in the free electron spectrum on the network. 
As soon as the magnetic field is present, k\ has to be chosen from [— 7r/a + 5 /a^~5 /a] U [+(5/a,7r/a — 5 /a] where 
cos((5a) = emoa;.(/)/v'l8, EmaxXf) denoting the largest eigenvalue of Ho at reduced flux /. Thus, the free electron 
spectrum exhibits an infinite number of non overlapping bands, for generic /. Gaps in the spectrum have two different 
origins. Most of them corresponds to change the internal motion of the particle along a single link, which amounts 
to turn n into n + 1. Other gaps are induced by the magnetic field, and already appears as gaps in the set of 
allowed A values (A spectrum). For the special value / = 1/2, the cage effect is manifested by the presence of a pure 

point-spectrum at energies £"„ = 2ma^ (±7r/2 ± a 4- 2?™)^, with a = arcsin(l/V3). Each of these levels is highly 
degenerate, reflecting all the possible cage states in the network. 

Since some experiments have been performed on superconducting wire networks, it is worth recalling briefly the 
connection between the A spectrum and the estimate of the superconducting critical temperature as a function of 
the external flux /. At the transition, the order parameter builds up from the nmpagating modes of the linear 
Ginzburg-Landau equation (written here in zero fleld for simplicity) along each wireQij : 

-|l + ^^ = 0- (31) 

The critical temperature is obtained from the largest value of ^ denoted by for which eigenvalues exist. ^* is 
known as the Ginzburg-Landau macroscopic coherence length and we have Tc{f) — Tc(0) cx l/(^*)^. To determine 
we follow the same procedure as for the free particle spectrum. We then have (fc*)^ = (^*)^ with k* G [0, 7r/a] and 
cos(fc*a) = Xmax.if)- ^or an infinite lattice, one has : XmaxXf) = £max.(/)/'\/l8 so that the Tc{f) curve is directly 
related to the edge of the tight-binding spectrum. This mapping has led to accurate comparisons between this simple 
model and the experimental datacl. However, for free boundaries system, such a simple correspondence no longer 
holds, since the coordination number can be different from 3 and 6 on the edges. For these reasons, our discussion of 
edge states will be given mostly in term of the linear problem ( |30| ) rather than in the tight-binding language. The 
main striking result is the appearance of very sharp edge states at / = 1/2. For those states, the amplitudes ipi^s are 
vanishing on most of the sites except on a finite width strip concentrated near the boundaries (see in Fig. [l^ ). 




FIG. 13. An example of edge states in a Ts wire network. We have put a black circle on each node for which the wave 
function has a non vanishing amplitude. 

For any shape of the boundary, these edge states are always dispersive (along the edges) . Note that in the tight- 
binding version, non dispersive edge states can also appear depending on the boundary shape. Another qualitative 
difference is that, in the tight-binding problem, the energy of these states always appears inside the main gaps of the 
infinite lattice spectrum whereas in the wire network, the A spectrum exhibits edge states with energy outside the 
bulk spectrum. We have displayed in Fig. ^ the spectrum of the finite lattice shown in Fig. |l^ and the infinite lattice 
spectrum. 
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1/3 1/2 2/3 
Reduced flux 



FIG. 14. Comparison of the A spectrum for the infinite lattice (right) and the finite size lattice displayed in Fig. |^ (left). 

The dispersive edge states are clearly visible. These have important consequences for superconducting networks 
since they offer the possibility to nucleate superconductivity on the edges of the sample, at temperatures slightly 
above the bulk T^iffiki. 

Notice that for specific boundaries, the dispersive edge states can run all around the sample as illustrated in Fig. |l^, 
providing possibly an interesting example of a non chiral (quasi) one-dimensional conductor. Some experiments on 
ballistic networks etched in a two-dimensional electron gas in GaAs riieterostructures have given some indications 
that edge states may have some influence in transport measurementE3. Indeed, the variation of the conductance 
through such a network with external field seems to depend on the geometry of the current pattern. On some 
samples, a stronger dip in the conductance at / = 1/2 has been observed when the current is injected directly in 
the bulk, in comparison to the usual setup where the current is injected and collected on the edges. In this latter 
configuration, a non vanishing conductance may be found for a perfect system thanks to these propagating edge 
statesllj. For illustration, we give an analytical expression for edge states (A spectrum) coresponding to the semi- 
infinite 7^ network displayed in Fig. |l^. For these specific edges, the only non vanishing amplitudes are A, C and D 
{B — 0). Therefore, we may choose a convenient gauge in the edge region as displayed in Fig. |l^. 















1 


1 >^ 


1 


1 





FIG. 15. A possible edge choice for the semi-infinite T3 network (same convention as in Fig. 

Assuming propagating waves in the edge direction with wave vector k, we must satisfy the following set of equations : 

3AA = C + 2Dcos(fca) (32) 



= 2Ccos(fca) - D 
3AC A 

3XD = 2Acos{ka) . 



(33) 
(34) 
(35) 



The corresponding A spectrum is given by A = ±i-\/l + 4cos^(fca). This example shows a dispersive edge spectrum 
A e [l/3,\/5/3] that lies below and above A = 1/V3 which is the maximum value of A for the infinite system. As 
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previously discussed, these states are thus more favourable than the bulk states for the nucleation of superconducting 
regions in a wire network upon cooling from high temperature. 

We have also investigated how these very narrow edg e states evolve as some perturbations away from the pure 
system at / = 1/2 are gradually introduced. On Fig. |6[ some edge states are represented for a pure system (a — 0) 
at / = 0.52 and for a disordered system (ct = 0.01) at / = 1/2. Both cases are very similar. The main feature is 
the spatial broadening of the wave function towards the center of the sample. However, we observe a rather smooth 
evolution of these states as the strength of the perturbation is increased. 



FIG. 16. An example of edge states for / — 0.52, cr = (left) and / — 1/2, a = 0.01 (right). The grey level on site i is 
proportional to the log 

Up to now, we have always consider that the electrons were independent. In the last section, we shall try to 
analyze the importance of the interactions between particles, in the lattice. Note that for simpler systems, the two 
interacting particle problem under magnetic field has already revealed some unusual interesting featureso. 



VI. THE INTERACTING CASE 



Of course, the many-body problem is certainly one of the most difficult to tackle. Here, our goal is to study, in 
a simple approach, the effect of the electron-electron interactions on the Aharonov-Bohm cages. We underline that 
this type of system provides a good starting point to understand the competition between localization and interaction 
even if the nature of the confinement is due to the magnetic field and not to disorder. In our case, the main departure 
from a realistic disordered model is indeed the preservation of the translation invariance in a localized system. In the 
following, we first focus on the two electron problem that already displays very interesting features, and we fix the 
reduced flux to its "critical value" / = 1/2. The system is described by the following Hubbard hamiltonian : 

Hnub. = ^ tij c\^^ Cj^a + m^i (36) 

(ij>,o-=Ta » 

= Ho + Hint. , (37) 

where c| ^ and Ci^^ denotes the creation and annihilation operator of a fermion with spin a respectively, rii^a- = c| o- '^i,'^ 
the density of spin a fermion on site i, and (. . .) stands for nearest neighbor pairs. Note that, since the particles 
considered here are fermions, the interaction term U is only efficient in the singlet sector where the orbital part of the 
wave function is symmetric. For simplicity, we completely neglect the coupling between the magnetic field and the 
electron spin. r— ■ 

Here, we shall mainly discuss the case of the 7^ lattice but we must mention the work presented in ref.ta on the 
chain of loops displayed in Fig. ^ For this system, we have shown that a Hubbard-like interaction term leads to a 
delocalization which is directly related to the emergence of dispersive bands in the two-particle spectrum. Of course, 
such a study is much more complex in the lattice, but, as we shall see, it is however possible to show that the same 
effect is present in this two-dimensional structure. 
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First, let us remark that if, at a given time, both electrons are far from each other, the local particle density is non 
vanishing in a cage that is simply the superposition of each individual cages. The interesting case thus arises when 
the electrons get close. We consider, for instance, an initial state where both electrons t and | are located on the 
same six-fold coordinated site, and compute the coefficients of the recursion chain for this two-particle wave function. 
With the site numbering given in Fig. ^, we denote this initial state by |(y5o) — |0, 0). More generally, a ket \i,j) 
represents a state where the particle with spin "f is located on site i and the particle with spin | is located on site 
j. Note that we must, in principle, work with the symmetrized ket but, as shown below, since our initial ket |0, 0) is 
already symmetric, it will automatically generate symmetrized kets. Applying the recursion algorithm in the gauge 



displayed in Fig. 17, one obtains : 

6 

HHub. \^o) = U\0, 0) + ^ \t, 0) + |0, z) (38) 

1=1 

1 ^ 

[/, 6i = 2V3, |^i) = ^^|z,0) + |0,z), (39) 
1 

73 

ai-0,52 = 2V3, 1(^2) = ^ ^ (41) 



ao 

1 ^ 

6 



u ' 



i?ff«fc.|'P2) = ^5]N,*)+62|^l) (42) 



i—i 



a2 



^3 = , 1^3) - I ^ E 1^' ^) - 6 E J^' ' (43) 



6 6 



Huub.Wz) = I ^E|i,i) + E|i,i-H6) - |i,i + 7) + K + 6,i> - + J (44) 



5 1 6 1 

a3 = ^, ^4 = 2^-, \^i) = —j=y2\i,i + (S)-\i,i + 7) + \i + Q,i)~\i + 7.,i). (45) 



Since [(^4) has a non vanishing amplitude on the second shell sites (i > 7) the two-particle wave packet can spread, 
so that the Aharonov-Bohm cages associated to |0,0) seems to be unbound or, at least, larger than the one obtained 
for |0) at C/ = 0. Actually, we have checked numerically that the 5i>5's are non vanishing so that the cage is really 
unbound. 
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FIG. 17. Aharanov-Bohm cage of a two-particle state initially localized on the central site, and its recursion chain (same 
convention as in Fig. For U 0, the black sites are reached. 

As for the chain of loops, the propagation is made possible by the existence of two-body bound states in the 
spectrum. The same analysis as in ref.L3 could be done in the lattice but the two-dimensional character of this 
later structure makes it more difficult. Nevertheless, if one focuses on the small interaction case {U <C t) it is possible 
to show the emergence of dispersive bands near the ground state energy. Indeed, one can, in this case, treat U as a 
perturbation on the infinitely degenerate level 2£_ without considering the higher energy levels. 

To solve this problem, we shall use the periodic gauge displayed in Fig. |l^ that provides real hopping terms. 




FIG. 18. A possible gauge choice for the Ts lattice at f — 1/2 (same convention as in Fig. g). 

Remark that the periodicity of the structure has doubled and we must now distinguish two different types of 6-fold 
coordinated sites denoted by (•) and (o). If one restricts the two-body problem analysis to the subspace corresponding 
to 2e_ (for t/ = 0), it is convenient to build the two-particle state basis as the tensor product of the one-particle 
orthogonal "cage" basis displayed in Fig. |l^. 
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FIG. 19. Non- normalized one-particle cage eigenstates associated to the two types of 6- fold coordinated sites (see Fig. [l8[ ). 
The amplitude on each site is given. 

Since the problem is invariant under a translation of the center of mass of the two particles, the basis of the singlet 
states sensitive to U can be expressed in terms of the following Bloch waves : 



l<Po(K)) = 

\vim = 

kuo(K)) = 

ku.(K)) = 

l^±v(K)) - 
k±w(K)) - 



r6A 



^ y e'K ' 



1 g.K., 

* re A 



' reA 
' rSA 



r + v,r + v) 



i',r + u)s 



r + v,r + V + u;5 



r, r ± v) 



r, r ± w) 



(46) 
(47) 
(48) 
(49) 
(50) 
(51) 



rSA 



where the origin O — (0, 0), u, v and w are represented in Fig. |l^ ; A is the lattice generated by u and v + w with 
periodic boundary conditions and containing Ns sites ; K is a wave vector lying in the first Brillouin zone associated 
to A. Note that the orthogonal basis used here to define the different vectors is different than the basis used in 
section [11. From now on, the ket |r, r')s represents the normalized space symmetric state in which one electron is 
in a cage state (eigenstate of the one-particle non interacting problem) localized around a six-fold coordinated site 
located in r and the other electron is in a cage state localized around a six- fold coordinated site located in r'. Because 
of the Pauli principle, these kets corresponds to antisymmetric spin wave functions and therefore to singlet states. 
Note that |r; r) is automatically symmetric so that we have omit the index S for this state. Naturally, if |r — r'| > ^/S, 
the electrons do not interact so that one may only to consider the 8 two-particle states defined above. Thus, in 
each irreducible representation indexed by K, one just has to diagonalize Huub. in the subspace generated by these 
vectors. Furthermore, as shown in the Appendix |A|, Huub. can be block diagonalized into two isospectral matrices. 
After shifting the energies so that for [/ = the eigenenergies are equal to zero (and not to 2e_), the eigenvalues (in 
U units) are given by the roots of the following characteristic polynomial : 



P(e,K) = e[10368e^ - 



288^) + e(258 + 2Q^i + 2u) + ^i- 3]. 



(52) 



/i = cos(K.u) -I- cos(K.v) -I- cos(K.w) (53) 
u = cos[K.(u + v)] + cos[K.(2u - v)] + cos[K.(2v - u)]. (54) 

Reminding that w = v — u, one can easily check that P is unchanged by the tranformation u ^ v, as it obviously 
should. 

It is readily seen in Eq.(|5^) that e = is always an eigenvalue for all K. The associated eigenvectors have a 
vanishing amplitude on the kets |(^o(K)) and |(pi(K)) that corresponds to states where both electrons are in the same 
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cage. The absence of dispersion for this family of states sug gest s that they may be generated by a set locahzed two 
particle singlets. As in the one particle problem (see Section III), it is even possible to exhibit such states where the 
center of mass is confined in a finite area. A possible choice is : 



\^Bound) - -= (|0, u)s - |0, -u)s - |0, v)s + |0, -V>5 + |0, w)5 + |0, -w)s) (55) 

The presence of these localized singlet bound states at zero energy seems to be a particular feature of this lattice. 

Apart from this non dispersive band, a close inspection of the polynomial P allows one to show that when K runs 
over an elementary cell of the reciprocal lattice, the roots of P describe two separate dispersive bands spreading from 
to 1/12 for the lower one and from 1/4 to 3/8. The corresponding eigenvectors are linear combinations of the 8 
vectors introduced above and are actually extended Bloch waves. These two-body bound states are delocalized and 
thus allow for a propagation of two-particle wave functions as we have already noticed using the recursion analysis 
(see if, Fig. The same perturbative approach in the degenerate ground state can be achieved for the chain of 
loopsEil. In this latter case, the spectrum, for U <^ t consists in two non dispersive bands with eigenenergies and 
1/16, and one dispersive band spreading also from 1/4 to 3/8. Here again these energies are measured in U units and 
shifted such that, for [/ = 0, the two-body ground state corresponds to a zero energy. We have displayed in Fig. ^ 
these two spectra. 

Two important differences arise between these two structures. First, there is no non dispersive singlet bound state 
at zero energy for the chain of loops, except the trivial ones involving two remote localized electron states. Non 
dispersive singlet bound states do exist but their energy is strictly positive and depends on U . Second, there is a gap 
between the U — Q ground state and the first dispersive band (bound state), by contrast to the Tij lattice. 
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FIG. 20. Perturbative low-energy two-particle spectrum for (a) the T3 lattice, (6) the chain of loops. 



The fact that all eigenenergies of the projected hamiltonian are positive is a simple consequence of the same 
property for the interacting part Hint. ■ Denoting by Q the projector on the lowest energy level of Hq, the lowest order 
degenerate perturbation theory proposed above amounts to diagonalize QHint.Q- For any two-particle eigenstate j?/'), 
we have {il}\Hint.\4') > if L/ > 0. In particular, this imposes all eigenvalues of QHint.Q to be positive (or vanishing). 
From this consideration, the appearance of non trivial (degenerate) zero energy bound states in the singlet sector is 
quite remarkable. 

Of course, it would be very interesting to characterize the properties of the system in the presence of a finite density 
of particles. In the following discussion, we shall always assume C/ <C so that we will only pay attention to the 
projected hamiltonian QHint.Q which is, as discussed above, a positive operator. Denoting by iVe the total number of 
electrons and by Ng the total number of sites, we introduce the density p = Ng/Ns that ranges from to 2 because of 
the spin degeneracy. In the non interacting case, for < p < 2/3, it is possible to put all the electrons, in the ground 
state of Hq. Thus, we will restrict our analysis to this interval. 

For p < 1/3, a possible choice for the ground state is to put each electron in a localized cage state (around 6- 
fold coordinated sites). In this case, the energy per particle is simply e_. At low density, if all occupied cages are 
further isolated from each other (no overlap), the interaction term Hint, does not affect these configurations which 
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are therefore ground states of QHnub.Q- These configurations have a huge spin degeneracy (2^=). One can also build 
ground state configurations containing some connected clusters of singly occupied cage states, provided these clusters 
are fully spin polarized. Obviously, the total spins of these clusters may be chosen independently from each other 
without any energy cost. Note that additional configurations with the same energy may also be obtained using the 
zero energy two particle localized singlet bound states discussed above. An interesting issue would be to determine 
whether this class of ground states exhausts all the possible ones, for p < 1/3. Note that if the density p lies beyond 
i/t/3 where vt is the site percolation density of the triangular lattice which is formed by the 6- fold coordinated sites, 
we have an infinite percolating cluster which leads to an infinite magnetization. 

As p approaches 1/3, the ground state thus becomes more and more polarized. If Ng denotes the total number of 
sites, the problem of finding the ground state of the projected hamiltonian becomes more interesting if 1/3 < p < 2/3. 
A simple class of eigenstates of QHu^b.Q is obtained if the magnetization is maximal. It corresponds to a total spin 
Smax. = — -/Ve/2. Indeed, because of the spin rotation invariance, such states can be built from states with 

N-f = Ns/3 and Ni = Ne/2 — Ns/3 after performing a global spin rotation. Here, and Ni denotes the total 
number of electrons with up and down spins respectively. The condition — Ns/3 means that all the single particle 
ground states of Hq are fully occupied by 1 spin electrons. For these states, the | spin electrons behave as completely 
localized free fermions with an individual excitation energy equal to U. To show this, we denote by ipr{i) — (j|r) the 
amplitude at site i of the state |r) which is the single particle cage eigenstates of Hq localized around the six- fold 
coordinated site located at site r, and we write : 

QH,nt.Q= X! "(l"3,r4|ri,r2) 4^^^4^ j^Cr2,iCri,T> (56) 
ri,r2,r3,r4 

with : 

u(r3,r4|ri,r2) = uY^ip*^{j)Lpl^{j)Lp^^{j)ipr,{j)- (57) 

3 

In expression (p6|), the fermions operator c\ ^ and Cr,cr create and destroy an electron with spin cr in a state |r). On 
a state where all the single particle ground states of Hq are fully occupied by up spins electrons, 4^ |Cri,| acts as 
(5r3,ri times the identity operator. It is then easy to check that : ^rU{YV4\r,r2) = U5r^^r2- As a consequence, for 
any state 1-0) such that 7V| = A^s/3, we have QHint.Q\4') = UNi\ip). At this stage, it not clear whether the ground 
state of the projected hamiltonian has the maximal value Smax. of the total spin for 1/3 < p < 2/3. Nevertheless, 
since this is true for p = 1/3 according to the previous paragraph, this generalization is plausible. The main question 
to be adressed is whether a J, spin may form a bound state with one or several "magnons" . By magnon, we refer 
here to particle-hole like excitations which destroy one t spin electron and create one J, spin electron. We leave this 
interesting but more complex problem for future investigations. 

To conclude this section, we may say that in spite of the presence of very low lying dispersive bound states for two 
particles, it is not so easy to turn the finite density system into a good conductor. However, we emphasize that some 
qualitative changes may arise for finite values of U, since a sizeable virtual occupancy of excited levels of Hq becomes 
then possible. 

VII. CONCLUSION AND PERSPECTIVES 

In this paper, we have studied three types of perturbations which affect the physical properties of a two-dimensional 
lattice embedded in a magnetic field that presents Aharonov-Bohm cages. Two of them, namely a weak Anderson-like 
disorder and the finite size effects do not drastically modify the main features of the ideal model. In particular, if the 
disorder is not too important, the single particle energy eigenstates remain strongly localized for half-integer fiuxes 
per elementary tile, and this independently of the disorder strength. This lasts until the disorder docs not introduce 
a strong mixing between distinct degenerate levels of the unperturbed tight-binding hamiltonian. By contrast, for 
generic values of the magnetic field, the amount of disorder-induced localization is very sensitive to the disorder 
strength. 

A more delicate situation is obtained in the presence of electron-electron interactions for the many-body system. 
We have shown how they partially destroy the single particle localization by the generation of extended two particle 
eigenstates. However, it is not clear at the present stage whether this mechanism is sufficient to induce a metallic 
behavior for finite electron densities and half a flux quantum per tile. This important question clearly deserves further 
investigations for the Hubbard model (studied here) and also for more realistic model of continuous narrow wires (few 
conducting channels), in connection to the experiments on two-dimensional electron gas. We expect to find partially 
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spin polarized ground states in a finite interval of electron density, and a fully polarized state when the number of 
electron per site equals 1/3. For slightly larger densities, an important question is to understand if the system create 
some spin textures resembling to skyrmions in quantum Hall effect for filling fractions close to odd integers. 

In this context, it would be very interesting to study interacting hard core bosons on this lattice, which would require 
numerical diagonalization on finite size clusters. Indeed, this model is closely related to the physics of Josephson 
junction arrays. In the limit of small capacitance islands quantum fluctuations of the order parameter phase are 
strong and the corresponding trend to localize Cooper pairs could be enhanced by the cage effect. In the other limit 
where the Josephson coupling dominates, it is not clear that the system at half-integer fluxes will develop a finite 
stiffness for spatial gradients of the order parameter phase. In this semi-classical regime, we expect a rather large 
ground state degeneracy as suggested by vortex decoration experiments on superconducting networks. 
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APPENDIX A: DIAGONALIZATION OF THE TWO-ELECTRON PROBLEM 



As explained in section VI, the non trivial part of the hamiltonian is given, in each irreducible representation 
indexed by K, by the subspace generated by the eight vectors : 

(l^o(K)), |^+v(K)), l^-v(K)) l^uo(K)), l^u, (K)), |^+w(K)), |^_w(K)), ^i(K))) . 

In this basis, the shifted hamiltonian _ff'(K) = (iJ/^ (K) — 2£_I^)/U (/§ denotes the 8x8 identity matrix) writes : 
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Note that since H'(K) is hermitian, we have only given its upper part. Now, let us introduce the following vectors : 
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that corresponds (up to a phase factor), to symmetric and antisymmetric combinations of the initials vectors. The 
hamiltonian does not connect the subspace generated by Bi = (|iy9(K))_ , |(y3v(K))_ |(/Pu(K))+, |iy9w(K))+) and the 
subspace generated by B2 = (|<^(K)) + , |(y9v(K)) + , |iy9u(K))„, |iy9w(K))_). In each subpace, H'(K) writes : 
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The characteristic polynomia P^i (e, K) and (^i K) associated to (K) and (K) respectively reads : 



P6,(e,K) = 



10368 

e 

10368 



[10368e^ - £^(3888 + 288^i) + £(258 + 20^1 + 2i/i) + a^i - 3] 
[10368£^ - £^(3888 + 288/13) + £(258 + 20/i2 + 21^2) + M2 - 3], 



with 



/ii = cos(K.u) + cos(K.v) + cos(K.w) 

^2 = cos(K.u) — cos(K.v) — cos(K.w) 

//I = cos[K.(u + v)] + cos[K.(2u - v)] + cos[K.(2v - u)] 

//2 - cos[K.(u + v)] - cos[K.(2u - v)] + cos[K.(2v - u)] . 



(A8) 
(A9) 



(AlO) 
(All) 
(A12) 
(A13) 



In the orthogonal basis chosen in section |l| where u = a (3/2,73/2), V = a(3/2,-\/3/2) and w = V — u, it is 
straightforward to show that i-'Bi(£, K) = Pb2{£,J^') if K = (fcijfe) and K' = {ki + 7r/3, A:2 — tt/v^). Thus, it is 
sufhcient to analyze the roots of the polynomial P{e, K) = Pgj^ {e, K) for all K belonging to an elementary cell of the 
reciprocal lattice. 
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